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We discuss the current status of our automatic perturbation theory program as applied to Fermilab Fermions. 
We give an overview of our methods, a discussion of tree level matching, and one loop results for the coefficients 
of the higher dimension kinetic operators. 



A dominant source of error in current lattice 
calculations is the errors due to finite lattice spac- 
ing a. One way to reduce these errors is to use 
improved actions. The general structure for an 
improved action is 



C = Cn 



. Cn{go,mo)Cn 



(1) 



The series here is an expansion in the dimension- 
ality of the various operators. Each of these terms 
comes with a new coupling constant c^. 

In order to use this type of action we must do 
two things, we must truncate the series at some 
specified order in n and we must calculate the new 
couplings somehow. Fixing the new couplings can 
be done in a number of ways, for example, the 
first step could be to make an expansion in pow- 
ers of m-o for light quarks, or powers of l/mg for 
heavy quarks. The remaining dependence on 
can be determined perturbativly. 

It is desirable to have a method of improvement 
which does not rely on light or heavy quark mass 
expansions. One such approach is the Fermilab 
approach |^ which orders the expansion opera- 
tors by dimension only (this amounts to a small 
momentum expansion). The unimproved fermi- 
lab action consists of dimension two and four op- 
erators. 

So ^ J dx'ip jmo + ^ ^ ^" -Pq - 2^" ^°^ 
+ CT-D-^Ajv' (2) 

For the definitions of the various derivatives we 
refer the reader to Q. This action includes 



an additional redundant dimension five operator, 
whose coefficient VsC can be tuned to remove the 
fermion doubling problem. 

To improve to 0{a^) a large number of of di- 
mension five and six operators must be added. 
The identification and tree level matching of these 
operators is discussed in |2|, we will not repeat 
it here. This report is concerned primarily with 
the one loop determination of the coefficients 
of the dimension six kinetic energy operators, 
{7 ■ D,!?^} and 7i-Df. For this matching we 
will set all of the other improvement coefficients 
(again see for a full list) to zero. This leaves 
the action, 

S^So + J dxt^[ci{j--D,D^} + c2-/^Df]^p. (3) 

For these kinetic operators the determination of 
the coefficients is very straightforward. For any 
lattice action, the quark energy can always be 
expanded in powers of the three momentum 
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By tuning the three couplings ci and C2, we can 
imposed various conditions on this action. The 
coupling C can be tuned to ensure that Mi = 
Al2, c\ and C2 can be tuned to ensure that A/2 = 
M4 and W4 = 0. The latter condition restores 
rotational symmetry to 0{a^). Because it only 
produces an overall shift in the zero of energy it is 
not necessary to tune M\= M2 using Rather, 
we'll just set ^ = 1 for convenience. 

With = 1 we have two remaining improve- 
ment conditions, M2 = and u'4 = 0. We 
impose these in perturbation theory. The tree 
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level calculation has been preformed in | 
ting Ts — C — ^, the tree level values are 
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In keeping with the Fermilab approach, these co- 
efficients have their full dependence on the bare 
mass. Notationally, bracketed superscripts on 
any quantity denote the order of the expansion 
in the bare coupling. For example, the rest mass 
Ml has the expansion 



Ml 



-9'aM['^ 



(7) 



The one loop calculations detailed below have 
been preformed using our automatic perturbation 
theory methods. The core of this method is the 
use of the Liischer - Weisz vertex generation al- 
gorithm , with some straightforward modifica- 
tions to include arbitrary quark actions. The ver- 
tex rules generated by this method can be used to 
construct Feynman diagrams, and VEGAS can be 
used to preform the loop sums. Additionally we 
use triple twisted periodic boundary conditions, 
with L — 200, as an infrared regulator. Further 
details on how these types of calculations are pre- 
formed can be found in Q . 

The automatically generated vertex rules can 
be used to check the tree level calculations in 
Q . Figure |l] shows the ratio of masses extracted 
from the automatically generated quark propaga- 
tor over a range of momenta. Correct tree level 
matching of cj"^ and C2°'' should give equation (^) 
with — Af4*''' over the full range of momenta. 
Clearly the improved action satisfies this, whereas 
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Figure 1. Tree level check on the energy, for a 
bare mass of one. The points show the ratio 



M, 



(0)/^^(0) 



over a range of momentum. The line 
labeled unimproved was generated with the bare 
action, the line labeled improved was generated 



with cf'' and C2'' equal to the tree level values, 
(I) and (I 

the unimproved action does not. We have also 
preformed this tree level checking for the tempo- 
ral one gluon vertex, verifying the expressions for 
the coefficients Ce, cg and cg in Q|. Work is cur- 
rently underway to fix and confirm the remaining 
tree level coefficients. 

We turn now to the one loop calculation of ci 
and C2. In order to fix these coefficients we must 
compute the quark self energy to one loop. De- 
tails of how to preform this calculation can be 
found in Our calculation mirrors |5|, apart 
from the different action with (H) and ^)), 
and the method of obtaining the Feynman rules. 
Expanding (jj) in powers of the bare coupling 
gives the one loop energy (recall Mj^"* — M^^^) 
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In terms of the quantities that appear here, the 
improvement conditions are, 



Mf =Mf, =0. 



(9) 
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Figure 2. One loop rest mass, over a range of 
input bare masses 

The one loop rest mass is easy to eval- 

uate numerically. Figure |^ shows its value over 
a wide range of bare masses. The limit of large 

(2) 

bare masses is approaching the value of = 
0.1681(4) reported in g. 

To compute ci and C2 we start by comput- 
ing using the tree level coefficients. Quan- 
tities computed with only the tree level values 
for ci and C2 will be denoted with bars. Setting 
Py = Pz = 0, and computing E'-^'> over a range 
of (small) Px allows and the combination 



w 
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to be extracted from fits to 
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the data. To separate AI^ from 0)]^' we com- 
pute i?^^) at fixed over a range ofpy. The cross 

term p^Py has coefficient ^,^^'^,„,\ which can be ex- 
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tracted from a fit to the numerical values of the 
one loop energy, as illustrated in Figure |3[ This 
procedure gives the following values (mg = 1), 



' = 0.08561(60) 
= 0.263(47) 



Mf ^ = 0.234(10) 

(2) 



-0.241(52) 



The coefficients c"i^ and c"2^ provide the coun- 
terterms that are necessary to correct for these 
one-loop contributions to the energy, in order to 
satisfy the improvement conditions (^ Hence we 
obtain the preliminary values (recall mo = 1): 



Figure 3. E^"^^ over a range of pi, the line is the 
best fit 

These values ensure the matching conditions are 

(2) (2) 

satisfied and give the value ' = Ml ' = 

0. 232.10). 

To conclude, this report illustrates the appli- 
cation of automatic perturbation theory to the 
Fermilab fermion action. These techniques can be 
used, either to do or to check, tree level matching, 
as well as one loop matching. We have presented 

f2) 

results for the one loop rest mass ' along with 
a first matching calculation for ci and C2. Work 
is now underway to extend the one loop matching 
presented here to more values of the bare mass at 
higher precision, and more coefficients. 
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